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Abstract 

Static vortices close together are studied for two different models in 2-dimen- 
sional Euclidean space. In a simple model for one complex field an expansion 
in the parameters describing the relative position of two vortices can be given 
in terms of trigonometric and exponential functions. The results are then 
compared to those of the Ginzburg-Landau theory of a superconductor in 
a magnetic field at the point between type-I and type-II superconductivity. 
For the angular dependence a similar pattern emerges in both models. The 
differences for the radial functions are studied up to third order. 
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1 Introduction 



Ever since 't Hooft |J and Polyakov found a monopole solution in the 
SU(2) Yang-Mills-Higgs theory, solitons in field theories have been studied 
extensively. Our understanding of monopole solutions has been greatly en- 
hanced by an existence proof for static solutions by Taubes [||] and the con- 
struction of monopole solutions started by Ward [|J. This process was not 
matched by quite the same progress in our understanding of the Abrikosov 
solutions of the Ginzburg-Landau theory, although one might have expected 
that the Abelian Higgs theory in 2+1 dimensions is actually simpler than the 
SU(2) Yang-Mills-Higgs theory in 3+1 dimensions. Again an existence proof 
was given by Taubes ||. However, only superimposed vortices can be de- 
scribed explicitly and no explicit construction of separated vortices is known. 
In this paper, we want to give the solution for two vortices close together in 
terms of an expansion in the parameters which describe the relative location. 

In Sections 2 and 3, we study a model for one complex field. Here the 
calculations are simpler than in the Ginzburg-Landau theory which is our 
second model. The first model has, however, some peculiar (unphysical) 
features. Assuming the most symmetric form in terms of angular dependence, 
only two smooth vortices can be superimposed, and when 'pulled apart', they 
develop a singularity at third order. In the Ginzburg-Landau model this does 
not happen. In fact, delicate cancellations take place to make the expansion 
smooth, at least up to third order. In this model the radial functions are 
given as solutions of certain linear ordinary differential equations. This is 
discussed in Section 4. 



2 Vortex solutions and zero modes in a simple model 

Our first model is a model 0|7]] for a pair of real fields 4> a (x), a,b = 1,2, or 
equivalently, for a complex field = 0i +i<p2- The Lagrangian density of the 
model reads 

C = d^dj^d^cP^ + (1 - |0| 2 ) 2 |0| 2 , (2.1) 
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where a,b = 1, 2 labels the components of the Higgs field and i,j = 1, 2 are 
the space indices. The square brackets mean antisymmetrization, 



d {i <f> a d j]( f> b = (d t r)(d^ b ) - (d j( p a )(d t ^). (2.2) 

We are working in 2-dimensional Euclidean space, i.e., the space indices 
can be raised and lowered without any change in the formulas. The indices 
which label the components of the Higgs field can also be raised and lowered 
without any change. In terms of the complex field <ft the Euler-Lagrange 
equation reads 

Wdji&Wty) = (1 - H 2 )H^(1 - I0DI0I- (2-3) 
Any solution of the equation 

2 det{ % ) = ±{l ~ 10,2)101 (2 ' 4) 



solves the equation of motion (2.3). Note that Eq. (2.4) is a first order 
equation whereas Eq. (2.3) is of second order. So we would expect that 
(2.4) is somewhat easier to solve than (2.3). For different types of models, 
this reduction of order was first introduced by Bogomolnyi ||. That is why 
we call Eq. (2.4) the Bogomolnyi equation here. Any solution of (2.4) also 
attains the lower bound in the following inequality, 

A = I Cd 2 x>^\Q\, (2.5) 
Jr 2 15 

where 

Q=~ [ te y (l -|0| 2 )|0|(9»(^0*)rf 2 x (2.6) 

is the winding number. Finally, all finite-action solutions actually solve the 
Bogomolnyi equation, so we do not miss out on any by concentrating on the 
first order equation. 

We now seek to attain a smooth finite- action solution of Eq. (2.4). For 

= f(r)e me (2.7) 
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Eq. (2.4) reduces to 



nf(r)f'(r) _ 1 



= -fi-F)f. (2-8) 

Since / — > as r — > (otherwise in (2.7) is not defined at the origin), we 
have 

„2 



/ = tanh— . (2.9) 
An v 1 

The solution in (2.7) with f(r) given by (2.9) is defined in the whole of 
R 2 and is clearly a C°° function in R 2 \ {0}. Since 

/»il-2exp — as r — > oo, (2.10) 
2n 

has the right asymptotic behaviour for a solution with winding number n. 
We still have to ensure that is C°° at the origin. There we use the Taylor 
expansion of /, 

oo o2fc/ 2 2fc _ D , r 2 r 2 I 2 

(2ife)! W ~4n 3 l 4n j + - 1 j 

where S fc is the k th Bernoulli number. We see that for n = 2 and only 
for n = 2, is a polynomial in x\ In this model, we have the (somewhat 
peculiar) situation that within the most natural ansatz (2.7) smooth finite 
action solutions exist only for n = 2, i.e., we only have a solution of the form 
(2.7) for 2 vortices. 

We have found the solution for two vortices sitting on top of each other, 
which we now denote by 0. To extend our study to two vortices slightly apart 
we consider = + 7, where 7 is very small, and we solve the Bogomolyni 
equation, linearized in 7. Equation (2.4) becomes 

2 d'y 2 , 2 d'y 1 

(/ cos6 l cos26 l H — /sin #sin2#)— '- + (f sine sin 26 + - / cos 6 cos 26) —— 
r ox r ox 1 

2 d'y 2 , 2 d'y 1 

-If sin # cos 20- -/cos#sin2#)— - (f cos sin 20 - -/sinflcos2#)— '- 

r ox r Ox 2 

= ^( 1 -3/ 2 ) (7 1 cos2fl + 7 2 sin2#) (2.12) 
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We find a 2-parameter family of zero modes, 



„2 



sinh — 

- r (r) = [a + P + i(a-P)]h(r) with h(r) = * . (2.13) 

cosh y 

These zero modes are C°° functions which vanish exponentially at infinity. 
By a rotation, one of the parameters could be removed and the vortices could 
be positioned on the x-axis, say. Since this does not simplify the calculations 
significantly, we will retain both parameters. Retaining the two parameters 
would also be necessary for a study of vortex scattering in the slow-motion 
approximation. This study is not done in this paper. 



3 The quadratic and cubic terms 



We now consider = + 7 + 5, and equate the second order terms in the 
Bogomolyni equation (2.4). This leads to the equation 

2 /M 2 /M 2 /M 1 rM 1 

-(/ cos 29^- + 2f sin 29^- - /' sin 29^- + 2f cos 29^-) 
r 09 Or o9 Or 

= {a 2 + (3 2 )fh\j 2 - 3) - ±fh 2 (3 + 1) [a 2 (cos 29 + sin 29) 2 

+2a/5(cos 2 29 - sin 2 29) + /3 2 (cos 29 - sin 29) 2 } 

+ (1 -3/ 2 )(5 1 cos2# + 5 2 sin2#) (3.1) 

with f{r) given in (2.9) and h(r) given in (2.13). 
With S of the form 

5 = a 2 F(r, 9) + 2a(3G(r, 9) + (3 2 H (r, 0), (3.2) 

we obtain the following equation for F(r, 9), 

9 (9F 2 (9F 2 

- (/' cos 26— + 2/ sin 20— - f sin 26— + 2f cos 20—) 

= ^ 2 (y - 3/) - y (3/ + j)(cos29 + sin2£) 2 

+ (1 -3/ 2 )(F 1 cos2# + F 2 sin2fl). (3.3) 
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To solve this equation we seek a solution of the form 

F = A(r)exp^-^ 2 (r)exp- 2e . (3.4) 

The ansatz (3.4) leads to two decoupled equations for fi and f 2 - In terms 
of the variable £ = r 2 /8, they read 

f + }(3f-l-f 1/^(1 -9f), (3.5) 

| + > (3/ »_ 1+ | )/l= _g (1+3 n (3.6) 

The general solutions to equation (3.5) is 

1 3 smh ^ smh ^ 
h = crf£ ( 2^sl?£ " c^sh£ + Cl) (3 ' 7) 

The function /i is a C°° function for < £ < oo. For £ — > 0, /i — > Ci holds. 
This implies that Ci = 0; otherwise F in (3.4) is not defined at the origin. 
Therefore, fi reads 

/i = 35inh / - (3.8) 
J 2 cosh 5 £ cosh 3 £ V ; 

The expansion of f\ near the origin is of the form 

oo oo 2 

k=i k=i ° 



Hence, the first term in (3.4) is a C°° function of x 1 and x 2 at the origin. We 
also see that f\ vanishes exponentially at infinity. So its contribution to 
does not change the winding number (2.6) which is a multiple of the action. 

A similar calculation yields a one parameter family of solutions to Eq. 
(3.6), namely 

f 2 = Sinh f - 3sinh3 / +C 2 ^. (3.10) 
2 cosh 3 £ 2 cosh 5 £ cosh 4 £ V ; 

In contrast to /i, all the solutions f 2 are acceptable. In fact, for all C 2 , f 2 is 
of the form 

oo oo 2 

f2 = Y,he = Y,b k (-) k (3.11) 

fc=l fe=l 5 
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near the origin, and therefore the second term in (3.4) is in C°°(R 2 ). The 
winding number and the action are also not altered because fa decays expo- 
nentially at infinity. 

The functions G and H in (3.2) can be found in the same way. If we put 
all results together, we obtain the second order terms, 

5 = (a 2 + (3 2 ) fa(r) exp l29 +t(a - t(3) 2 fa(r) exp^ 29 , (3.12) 

where fa and fa are given by (3.8) and (3.10), respectively. 

To find the cubic terms, we consider = + 7 + 5 + e, with 7 given in 
(2.13) and 5 given by (3.12). We set (3 = and concentrate on 

e = a 3 I(r,9). (3.13) 

For the Bogomolnyi equation to hold, I must satisfy 

2 dT 2 dT 2 dT 1 8T l 

- (/' cos 29— + 2f sin 29— - f sin 29— + 2f cos 29—) 
r 09 Or 09 or 

+h' (2 fa cos 29 + 2 fa sin 29) + h' (2 fa sin 29 + 2 fa cos 29) 
= -Sfil 1 cos 29 + I 2 sin 29) - fa(cos 29 + sin 29) 
-3//i(cos 20 + sin 20) (A - 2 fa cos 20 sin 20) - 3(cos 29 + sin 29)h 2 

+(J 1 cos 20 + J 2 sin 20) + - [fa (cos 20 + sin 20) - fa (cos 20 + sin 20)] 

2 

h 3 h 3 
+— (cos 20 + sin 20) 3 - —(cos 20 + sin 20) 
2 / 2 

- i (cos 20 + sin 20) (fa -2 fa cos 20 sin 20) + ^— (I 1 cos 20 + 1 2 sin 20) 3 (3.14) 
/ 2 J 

To solve equation (3.14) we seek a solution of the form 

/ 1 = ^i(O+^(O(cos40-sin40), 

I 2 = 9i(0 ~ ^ 2 (O(cos40 + sin 40). (3.15) 
This implies that g\ and g 2 must satisfy the equations 

df + (V " = ~ ^ " 6 ^ + " 2f - AP ~ 4/ ' (3 - 16) 



d£ 7 " J ' fdf^ 2/2 
The general solution to Eq. (3.17) is 



2 4/3 4 /- 



sinh£ 5sinh 3 £ / sinh 2 £ 3sinh 4 A sinh 3 £ 

9^ = ? t h Co a TT~ + L-3 f . 3.18 

y 4 cosh 5 £ 4 cosh 7 £ V 2cosh £ 2 cosh 6 £7 cosh 5 ^ 
All solutions (3.18) decay exponentially at infinity. For r — > 0, however, 

S 2 (r) = ^r 2 + ... (3.19) 

Hence, / in (3.15) is not a C°° function on R 2 . Our expansion gets singular 
at third order for the ansatz (3.15). In the next section we will discuss a 
realistic model in which a similar pattern emerges but no singularities occur. 



4 Abrikosov vortices 



The Ginzburg-Landau theory of a superconductor in a magnetic field in di- 
rection z is given by the Lagrangian density 

c = + ^(A0)p>r + ^(i r - 1) 2 , (4.i) 

where is the complex Higgs field, and D^cj) = di4> — iAi(p and = diAjdjAi 
in terms of the gauge potentials A iy % = 1,2. The Euler-Lagrange equations 
are 

D t D^ = ±<j>(l- | <P | 2 ), = l -[<P(Di<PY - <p*D^] (4.2) 

In the special case A = 1 it can be shown || that all finite action solutions 
of Eq. (4.2) satisfy the first-order Bogmolnyi equations |5J, 

F 12 = ^(1- |0| 2 ), D l( f> = -iD 2 <p. (4.3) 



It has also been shown || that a 2n-parameter family of solution of (4.3) 
exists with winding number 

n = ±- / F 12 d 2 x. (4.4) 

This family describes n vortices sitting at n position in space. 

Even for n vortices sitting on top of each other, the solution is not known 
explicitly in terms of elementary functions. It is known |i~0 |, however, that 
this solution is of the form 

4> = f(r)e me , Al = -^e^, (4.5) 

where / and a satisfy 

rf - n(l - a)f = (2n/r)a + f - 1 = (4.6) 

and 

/(0) = o(0) = 0, fim f(r) = fim a{r) = 1. (4.7) 

In the following, we restrict our attention to n = 2 and use the solution 
(4.5) as the zero order term in an expansion in the separation parameters. 
The first order terms are given by the two zero modes describing the sepa- 
ration of the votices. These were found by Weinberg [OH. Using his results 
we can write, up to quadratic terms, 

= f e 2ie + 2(a + i/3)kf + a 2 ip + a(3(f) + (3 2 X + ■ ■ ■ , (4.8) 

A x + iA 2 = iye l6 - 2i{a + i(3)(k' + y)e~ ld 

+ a 2 (5i + iB 2 ) + af3(d + iC 2 ) + p 2 (D x + iD 2 ) + ... (4.9) 
Here the radial function k(r) satisfies 

k" + h' - (f 2 + ^)k = 0, (4.10) 

with 

\mvr 2 k = 1, lim Mr) = 0. (4.11) 
Our task is to determine ip, <p, x, Bi, Ci, Di, which are functions of r and 9. 
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Equating the a 2 -terms in the Bogomolnyi equations (4.3), we obtain 

2n 2k 
(<9i + id 2 )4> + —i)e s - tf(B 1 + iB 2 )e 2td = 4kf(k' + y)e~ ie , (4.12) 

d 1 B 2 - d 2 B 1 + \{Me~™ + f^e 2ld ) = -2k 2 f. (4.13) 

A Fourier expansion with the minimal number of nonzero terms leads to the 
ansatz 

iP = g(r)f(r)e™ + ~g(r)e-™, 

B l + %B 2 = b(r)e ld + tb(r)f(r)e- 3ld , (4.14) 

and to equations for g(r),g(r), b(r) and 6(r). The equations for g{r) and b(r) 
read 

= I±^ 6 _ 6 ' 6 = _ lh '_ ( 415 ) 

r 

The functions g(r) and 6(r) must satisfy the equations 

9" + V - A = 2A; 2 / 2 , (4.16) 
r 



„» + i 6 < _ (i±Z! + l + 4» + 4«' = , + 2_k 

r 2 r l r 

Equating the a/3-terms and the (3 2 -terms in the Bogonolnyi equation (4.3), 
we obtain equations for and Cj, and for x and Dj respectively. These 
equations, which are very similar to equations (4.12) and (4.13), can again 
be solved by functions with the same ^-dependence as in (4.14) but with 
slightly different radial functions. Collecting all results, we can write the 
solution, up to quadratic terms, in the form 

= fe 2ie + 2(a + ifi)kf 

+ (a 2 + (3 2 )gfe™ + (a + i[3) 2 (±±^b - b')e~™ + ... 

r 

A x + iA 2 = tye 19 - 2i(a + if3){k' + y) e - ld 
- i[a 2 + P 2 )g'e ld + i(a + i^fbje^ + ... (4.18) 
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It remains to be shown that the quadratic terms in (4.18) are C°° func- 
tions on R 2 which do not change the action (and the winding number). To 
this end we use the asymptotic expansions of /, a and k at zero [|T2|], 

f{r) = f 1 r 2 + l -f l r A +... ) a(r) = ir 2 -l/ 2 r 6 + . . . , k(r) = r~ 2 +k ir 2 +. . . , 

(4.19) 

where f\ = .236 and k\ = —.025 from the numerical analysis. We find that 
the solutions of (4.16) and (4.17) have the following expansions at the origin, 

g(r) = g-i logr +g t + ^f 2 r 2 + ... 

b( r ) = b_ ir - 1 + b x r + (-61 - 2/iA;i)r 3 + . . . (4.20) 

8 

The higher order terms in g(r) are even powers of r, whereas the higher order 
term in b(r) are odd powers of r. Hence, the quadratic terms in (4.18) are 
C°° near the origin if and only if = b-\ = 0. So far the constants g\ and 
b\ are arbitrary. 

For large r the functions /, a, and k have the following asymptotic be- 



havior 1 2 



/(r) = l + A(r)e- r + ..., 

a(r) = l + a!(r)e- r + ..., (4.21) 

k(r) = ki{r)e~ r + . . . , 

with coefficient functions which are polynomially bounded. This leads to the 
existence of exponentially decaying solutions which asymptotically are of the 
form 

9(r) = gi(r)e- r + ..., b(r) = ~b x {r)e- r + ... (4.22) 

Here g\ and b\ are polynomially bounded. 

By numerical integration, the coefficients g\ and b\ which lead to an 
exponential fall-off at infinity, are found to be g\ = —.144 and b\ = —.026. 
The existence of such functions can be explained analytically as follows: 
Equation (4.16) shows that for positive g x , g cannot have a maximum for 
any r. So the function diverges exponentially. For very small g x , the term 
on the right-hand side of (4.16) will force the function to cross the r-axis, 
and then, as before, diverge exponentially. For very large negative g x , the 
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third term in (4.16) will force g to go through a maximum for large r. After 
that, the function cannot have a minimum and must go to minus infinity. 
Because of the continuous dependence on the initial data, we have an open 
set of data for which g crosses the r-axis, and an open set of data for which 
g goes through a maximum below the r-axis. Therefore, we have at least one 
value of <?i for which the function does neither. This function must converge 
and does so to zero, exponentially. 

A similar argument explains the existence of an acceptable solution b(r) 
to Eq. (4.17). The right-hand side of that equation is positive. So again 
b cannot have a maximum above the r-axis. Also, for very small negative 
b\, the right-hand side will force b to go through a minimum and then cross 
the r-axis. For very large negative bi, the third term in (4.17) prevents b 
from going through a minimum. In between these two possibilities we find 
the desired solution which goes through a minimum but does not cross the 
r-axis. Such a solution must decay exponentially. 

The cubic terms can be calculated in the same manner. We find, at third 
order, 

= . . . + (a + t(3)(a 2 + P 2 )fh + (a + t[3) 3 (-c> + ^±^c)e~ Aie + ..., 

A x +iA 2 = ... 
2 2k 

+ i(a+iP)(a 2 + p 2 )[-ti h + 2g(k' + — )+2kg']e~ %e + i{a + if3f fce~™ + . . . 

r r 

(4.23) 

The new radial functions, h(r) and c(r), satisfy the equations, 

h" + -ti - (f + \)h = Ak'g' + 2fk(2fk 2 + 3fg + ^^b - b'), (4.24) 



c - + i c --(l±/! + 9 + 12 ° + 4 °V = M/V . 2( v + »)(i±!V t 0. (4.25) 
t 2 r r r 

Near the origin, Eq. (4.25) has a series solution in powers of r 2 of the 
form 

h(r) = f 2 + hr 2 + h 2 r 4 + ... (4.26) 

The constant term is given in terms of the coefficient fi of the leading term in 
the expansion (4.19) of f(r). The form of this term leads to the cancellation 
of the r _1 -terms in the radial function multiplying e~ td in (4.24), and thus 
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ensures that this term in (4.23) is C°° on R 2 . The series in odd powers of r 
for c(r) which solves Eq. (4.25) near the origin, is 

c(r) = cir 3 + c 2 r 5 + . . . (4.27) 

The form of the series solutions at the origin guarantees that the cubic terms 
in (4.23) are C°° functions on R 2 . For large r, Eqs (4.24) and (4.25) have 
exponentially decaying solutions. 



5 Conclusions 



Our expansions show a simple ^-dependence in terms of trigonometric func- 
tions. In both models, the expansion of exhibits the following pattern: 



e 2ld 



e 0i8 



e -2i0 e 2i9 
e -4i6 e 0t6 
3 -6«6» e -2i0 e 2i0 

e -4i6 e 0t6 



Here the first line gives the 6 dependence of the zero order term; the second 
line gives the first order term, and so on. We get a similar triangular pattern 
for the 8 dependence of A\ + 1A2 at any order. For the radial functions we 
find differences between the two models. In the model for one complex field, 
the radial functions can be given explicitly in terms of exponential functions. 
However, for the angular dependence (3.15), a singularity occurs at the origin. 
(We have found no solution to (3.14) which is not of the form (3.15); we have 
no proof that there is none.) 

For the Ginzburg-Landau theory on the other hand, the expansion is 
smooth, at least up to the order to which we carried out our calculations. In 
this model the radial functions are not given in terms of well-known functions. 
Having used the technique to calculate the terms up to third order, it is quite 
clear how to proceed to any order, and also how to proceed in the case of 
more than two vortices. We expect these expansions to converge for small 
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separation parameters in the physical Ginzburg-Landau model. However, we 
do not have an estimate of the radius of convergence. 
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